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ABSTRACT 

We estimate the total dust input from the cool evolved stars in the Small Magellanic Cloud (SMC), 
using the 8-^m excess emission as a proxy for the dust-production rate. We find that Asymptotic Giant 
Branch (AGB) and red supergiant (RSG) stars produce (8.6-9.5) x 10~ 7 M yr _1 of dust, depending 
on the fraction of far-infrared sources that belong to the evolved star population (with 10%-50% 
uncertainty in individual dust-production rates). RSGs contribute the least (<4%), while carbon- rich 
AGB stars (especially the so-called "extreme" AGB stars) account for 87%-89% of the total dust 
input from cool evolved stars. We also estimate the dust input from hot stars and supernovae (SNe), 
and find that if SNe produce 10~ 3 Mq of dust each, then the total SN dust input and AGB input 
are roughly equivalent. We consider several scenarios of SNe dust production and destruction and 
find that the interstellar medium (ISM) dust can be accounted for solely by stellar sources if all SNe 
produce dust in the quantities seen around the dustiest examples and if most SNe explode in dense 
regions where much of the ISM dust is shielded from the shocks. We find that AGB stars contribute 
only 2.1% of the ISM dust. Without a net positive contribution from SNe to the dust budget, this 
suggests that dust must grow in the ISM or be formed by another unknown mechanism. 
Subject headings: stars: AGB - ISM: dust, extinction - galaxies: Magellanic Clouds - stars: supernovae 



1. INTRODUCTION 

Dust in galaxies plays an important role in allowing 
molecular clouds to cool sufficiently to form stars and 
planets. Dust is known to form in the atmospheres of 
stars or in explosive/eruptive events such as novae and 
supernovae (SNe). There is no known mechanism for 
forming dust in the interstellar medium (ISM), though 
it may be possibl e to grow dust in the ISM from existing 
dust grains (e.g.. |Pwe5[l998t IDraind l2009L The rate of 
grain growth within the ISM can be estimated by know- 
ing the dust-injection rate from stellar sources and the 
dust lifetime within a particular galaxy. 

Asymptotic Giant Branch (AGB) stars are impor- 
tant dust creators in galaxies. Rece nt work, how ever, 
has questioned the long-held belief ([Gehrzl I1989D that 
they are the primary dust factories. For the case of 
Jl 148 16+525150, one of t he most distant quasars (z = 
6.4), IValiante et all ((2009) argue that AGB stars must 
be the source of the obse rved dust (2 x 1O 8 M ), but 
iDwek fc Cherchnefil (pOll argue that SNe can produce 
the observed dust if the star formation history is assumed 
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to have long periods of very low star formation rates. SNe 
are generally thought incapable of producing enough dust 
to ac count for the amoun t seen in the ISM (<0.01 Mq 
each; lAndrews et aL1l2011l, and references therein). Sur- 
prisingly. iMatsuura et al.l(|2011l ) recently detected nearly 
1 Mq of dust around SN1987A (see also lLakicevic et all 
|2011[); the equivalent dust- injection rate requires 10 7 - 
10 10 dusty AGB stars. However, it is unknown whether 
SNe have a net positive or negative impact on dust pro- 
duction, as they are also efficient dust destroyers. 

The only way to compare the total dust input from the 
evolved stars in a galaxy to the total dust budget is to de- 
tect the entire population of dusty stars at infrared (IR) 
wavelengths and estimate the dust-injection rate of each. 
These sorts of global measurements are difficult in our 
Galaxy owing to obscuration by the Galactic Plane and 
in Local Group dwarf galaxies owing to li mited sensitiv- 
ity and resolution (cf. iBover et al.ll20 09b). Only for the 
Magellanic Clouds has the appropriate wavelength cov- 
erage, sensitivity, spatial coverage, and resolution been 
achieved to attempt to derive the total dust budget of 
a star- f orming galaxy. I n the Large Magellanic Cloud 
(LMCh lMatsuura et al.l (|2009D find that AGB stars and 
SNe combined account for only 3% of the ISM dust, 
though this could be much higher if other SNe remnants 
have large dust reservoirs like that around SN 1987A. 

Here, we estimate the current global dust production 
from the entire population of AGB and red supergiant 
(RSG) stars in the Small Magellanic Cloud (SMC). With 
a me tallicity that is 2.5-3 times lower than the LMC 
(e.g.. lRussell fc Donital[l992T: ILuck et al.lll998T >. the SMC 
is more representative of high-redshift galaxies. Com- 
bined with estimates of dust production in SNe, we esti- 
mate the total dust budget of the SMC and compare the 
results to measurements of the ISM dust mass. 
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1.1. Estimating Dust- Production Rates in Evolved Stars 

The effect of metallicity on the amount and type of 
dust species produced by an evolved star is not yet well 
unde r stood. Previous stud i es (e.g.. IGroenewegen et al.l 
120071: iMatsuura et al.ll2007t ISloan et al.ll2008D find that 
more metal-rich stars do appear to produce more oxygen- 
rich dust. Presumably, this is because O-rich dust pro- 
duction is limited by the metallicity-dependent availabil- 
ity of oxygen and silicon, while C-rich dust production 
is limited by the mostly metallicity-independent produc- 
tion of carbon by the star itself. Metal-poor stars are 
thus more likely to produce C-rich dust, as any oxygen is 
quickly tied up in CO molecules after dredge-up, leaving 
an excess of carbon. However, there is evidence that even 
carbon stars p r oduce less dust at lo w metallicity (e.g., 
Ivan Loonll200Ct Ivan Loon et al.ll2008f ) . possibly due to a 
lack of nucleation seeds. 

To accurately estimate the total amount of dust 
production around a star, detailed radiative transfer 
modeling of a well-sampled spectral energy distribu- 
tion (SED) is necessary. In addition, the inclusion 
of an IR spectrum to determine the dust species is 
ideal. Studies using this method have been carried 
out for subsets of the dust-producing population in 
several galaxies and clusters spanning a wide range of 
meta lli city (e.g., Ivan Loon et al.l 120051 : IMatsuura et al.l 



20071: IGroenewegen etaTj | 2007t iBover et al 
Groenewegen et al.l I2009alrbl : ILagadec et al 



McDonald et al.l 120091 l2011alf . However, it is currently 
unfeasible to use this method to measure the global dust 
input of an entire population of dust-producing stars 
since detailed radiative transfer modeling of thousands 
of individual stars is a prohibitively long process. 
One must thus depend on photometric techniques for 
estimating the dust input, usually in the form of an 
IR color analysis since IR colors generally scale with 
the mass-loss rate (e.g.. | Groenewegenll2006t ISloan et al.l 
l200l IMatsuura et al.1 120091 ) . In this work, a full SED 
from optical to IR wavelengths is available for each 
SMC evolved star, allowing us to use the more physical 
approach of computing the IR excess over the stellar 
photosphere as a pro xy for the dust production (cf. 
Srinivasa n et al.1 I2009D rather than relying on the IR 
colors alone. 



2. DATA & ANALYSIS 

The photometric data used in this study are from the 
Spitzer Legacy program "Surveying the Agents of Galaxy 
Evolution in the SMC", or SAGE-SMC. Images were ob- 
tained at 3.6, 4.5, 5.8, 8.0, 24, 70, and 160 fim covering 30 
deg 2 , including the Tail, Wing, and Bar of the SMC. Op- 
tical to near-IR photometry fr om the Magellanic C louds 
Photometric Survey (MCPS: iZaritskv et al.1 12002D . In- 
fraRed Survey Facility (IRSF; iKato et al.ll2007l). and the 
2-Micron All Sky Survey (2MASS: ISkrutskie et al.l[2006h 
were matched to the Spitzer photometry, and the full cat- 
alog is available for download through the Spitzer Science 
Center 8 . For details on the conte nt of the cata l og an d 
the nature of the observations, see IGordon et all (|20H . 

Throug hout, we adopt Ay = 0.12 mag and E(B — V) = 
0.04 mag (jSchlegel et al.lll998l: lHarris fc Zarltskvl l200l 



to account for interstellar reddening. We use the extinc- 
tion law from iGlasi (|1999l ) for o ptical to near-IR bands . 
Extinction in the Spitzer bands (jlndebetouw et al]|2005l) 
is negligible (A x = (6-8) x 10" 3 mag, for A = 3.6- 
8 Atm). We assume the distance to the SMC is 60 kpc 
(|Cioni et al.ll2000l: Keller fc Woodll2006l) . 

AGB and RSG star candidates were sel ected photo- 
metric ally from the catalog, as described by IBover et al.1 
(|2011h . AGB candidates were separated into carbon-rich 
(C-AGB), oxygen-rich (O-AGB), extreme (x-AGB), and 
anomalous O-rich (aO-AGB) sour ces. The x-AGB sam- 
ple is dominate d by carbon stars (Ivan Loon et al.l 119971: 
van Loonl 120061 : Ivan Loon et~aTl 120081 : IMatsuura et al.1 



20091 ) . though it likely includes a small number of ex- 
treme O-rich sources. The aO-AGB sources are a sub- 
class of the O-AGB candidates with redder J — [8] colors 
than the bulk O-AGB popu lation at the same magni- 
tude. See IBover et al.1 ()201lD for a detailed description 
of their IR properties. Since the nature of th e aO-A GB 
sources is unknown and since ISrinivasan et al.1 ()2009f ) did 
not distinguish between the aO- and O-AGB sources for 
the LMC analysis, we lump them together for the mass- 
loss analysis (Sect. |2~!3|) . 

Some cross-contamination between each stellar type 
and from young stellar objects (YSOs) and other in- 
terlopers is likely, but we expect it to be minimal (see 
IBover et al.l 120111 ). Evolved star candidates showing 
F24fim > fs/im are more likely to include contamination 
fr om YSOs or planet ary nebulae (Far-IR - FIR - objects 
in iBoyeret al. 2011), so these are considered separately 
from the AGB and RSG samples. 

2.1. IR excess 



http: //data, spitzer . caltech.edu / popular / sage-smc / 



We followed ISrinivasan et al.l ()2009l ) to estimate the 
8-fj.m excesses (Y 8 _ pm ) of the AGB and RSG stars. 
The 24- fim excess was not considered here because 
<20% of the O-rich st ars are detected at 24 /xm. Sec 
ISrinivasan et al.l (|2009f ) for a comparison between the 8- 
and 24-/zm excesses in the LMC. To compute A 8 _ Mm , we 
first fit the photospheric emission at optical and near-IR 
waveleng ths to plane-parallel C-rich COMARCS mod- 
els from IGautschv-Loidl et al.1 (|200l for the C-AGB 
stars and the sp h erical O-rich PHOENIX models from 
IHauschildt eFatl (fl999l) for the O-AGB, aO-AGB, and 
RSG stars. We chose one model for each type of star 
that best fit SEDs with little or no dust. The best- fit 
model photosphere was then scaled to the f/-band flux 
to estimate the IR excess at 8.0 /xm. Luminosities were 
determined by a simple trapezoidal integration of the U- 
band to 24- /im flux. 

The x-AGB stars and FIR objects are so heavily ob- 
scured in the optical as to make it impossible to fit 
the stellar photosphere. For these sources, we there- 
fore assumed that the flux in the IR is completely dom- 
inated by the IR excess. This is a reasonable assump- 
tion since we expect <15% of the mid-IR flux to come 
from the stellar photosphere itself in x-AGB stars (for 
J- [3.6] > 3.7 mag). 

The dust mass-loss rate, or dust-production rate 
(DPR), is expected to scale with luminosity as tL (e.g., 
Ilvezic fc Elitzurlll995l ). We therefore expect the DPR 
(D) to increase with luminosity (hence evolution), pro- 
vided optical depth (r) does not decrease with evolution 
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faster than L^ 1 . Figures 7, 8 and 9 fromlSrinivasa n~et al.l 
( 2009) show that the IR excess also scales with luminos- 
ity for AGB stars, making the excess a good proxy for the 
DPR. The relationship between luminosity and ^s-^im in 
the LMC agrees with what we find in the SMC within 
the uncertainties. 

The 8-fj.m excess may not scale well with the DPR for 
the dustiest O-rich stars that show silicates in absorp- 
tion at 10 /im rather than in emission. We expect these 
to be rare in the SMC, where the dust fraction in O-rich 
stars is low owing to its lower metallicity, and the smaller 
population results in fewer stars in this short-lived evo- 
lutionary phase. 

To avoid false detections, iSrinivasan et al.l (|2009l ) de- 
fine a threshold for reliable excesses in terms of data qual- 
ity - only excesses with relative uncertainties less than 
1/3 (hereafter, ">3-u excess") were used in the analysis. 
We follow this convention in the current paper. 

2.2. Derivation of Dust- Production Rates 
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Fig. 1. — Color-magnitude diagrams (CMD) of evolved stars fit 
with the GRAMS models. The underlying Hess diagram (gray) 
represents the full SMC catalog. Upper Panel: Near-IR CMD. 
Solid lines show the approximate d ivision between su pergiants, O- 
rich sources, and C-rich sources (cf. Boycr ct al. 2011). The dotted 
line marks the tip of the Red Giant Branch (TRGB). Lower Panel: 
Mid-IR CMD. Stars were selected for fitting if they have >3-<r 
excess and either an IRS spectrum or AKARI photometry exists 
to better constrain the SED. The fitted stars are a representative 
sample of the SMC evolved stars. Example SEDs are shown in 
Figure [2] 

In order to convert the IR excess to the DPR, we 
require a set of ^taxs with known DPRs in the SMC. 
iGroenewegen et aL (2009b) performed detailed radiative 



transfer modeling of several stars in the SMC, but these 
include mainly x-AGB stars and the dustiest O-AGB 
stars. We therefore select a subset of each type of evolved 
star and model their SEDs to determine their DPRs. We 
select sources: 

1. that represent SMC AGB and RSG stars over the 
full IR color space (Fig. [TJ , 

2. that have an 8-/iui excess with quality >3c, 

3. whos e photometr i c clas sification as O- or C-rich 
from iBover et al.1 ()201lD does not contradict the 
spectr osc opic classification from [v an Loon et al.l 
(|2008| ) or IGroenewegen et all (|2009bl ). and 

4. that have additional AKARI photometry and/or 
Spitzer InfraRed Spectrograph (IRS) spectra to 
create a better-constrained SED. IRS spectra are 
from the SMC-Spec Sp i tzer p rogram (P.I. G. 
Sloan); see iKemper et al.l (|2010l ) for a description 
of IRS data red uction. AKARI photometry is from 
llta et all IpOloh . 

To find the DPRs of this subset of evolved stars, 
we use the Grid o f RSG and AGB Mod elS (GRAMS; 
iSargent et aI1l2011t ISrinivasan et al.ll2011l) . These mod- 
els cover the full range of stellar and dust properties rel- 
evant to RSG and AGB stars and were developed to 
reproduce the IR c olors of evolved sta r s in th e LMC. 
GRAMS fits to the IGroenewegen et alj (|2009bT ) carbon 
stars produce DPRs that are systematically lower by a 
factor of 2-4 due to a difference in opacities of the amor- 
ph ous carbon dust used in the two studies (see Fig. 11 
in ISrinivasan et alJl20ll . The GR AMS O-rich fits do 
not sho w a systematic offset from the IGroenewegen et al.l 
(2009b) values; in fact t here is an overall agr eement be- 
tween DPRs (Fig. 14 in ISargent et al.|[2^Tll ). However, 
it is still possible that fits for individual stars differ in 
DPRs by a factor of up to 6, especially at lower DPRs 
(this discrepancy is also due to different optical constants 
for silicate dust). See Section l2~5l for a discussion of un- 
certainties. 

We fit the SEDs of the selected SMC sources with 
GRAMS models using a simple chi-squared routine. The 
chi-squared calculation includes AKARI photometry, 
where available. IRS spectra are not included in the 
chi-squared computation, but are used to confirm by eye 
whether the best-fit model is a good match. Sources with 
poor GRAMS fits are excluded from the sample, leaving 
12 O-AGB, 16 C-AGB, 13 x-AGB, 12 aO-AGB, and 14 
RSG stars. We show a sample of these stars in Figure [2] 
and their DPRs are listed in Table [1] 

GRAMS assumes a wind expansion velocity of 
10 km s _1 for all stars. Assuming that t> cxp is 10 km s _1 
for a star with L = 30 000L Q in the LMC, we scaled 
the DPR from GRAMS according to the following: D oc 



L QJ 



^cxp 



and 



dust ratio 



(cf. Ivar 



0.5 



van Loon et al.l 



where tp is the gas-to- 
2006). We assume that 



%j) sca les with metallicity (jvan Loonl l2000: Marshal l et al 
' 2004 1 ) Zrmo = 0-2 \ Z Q , and ->p & = 200 (e.g- lKnapp et al 
1993t lKnapp| [2001). so that ip SM c = 1000. However, we 



note that the gas-to-dust ratio metallicity dependence re- 
mains highly uncertain, and may not be the same for O- 
rich and C-rich sources. It has been suggested that C-rich 



Boyer et al. 



J004326.48-732643.3 



J01 0029.38-722257.6 



TABLE 1 

Dust mass-loss rates for stars fit with GRAMS models 



Spitzer Designation 3 


Name 


Type 


log(D) 
(M /yr) 


TOnO/TO T A ^ A 

J003656. 74— 722517.4 


MbA blVLO 091 


x-AGB 


—9.0 


Ton 1 tr o 7o 71 ATnn o 

J004o50. 78 — 714739.2 


MbA. bMO 200 


x-AGB 


—8.8 


TAA J CI f r\ AT TO O P" O O T 

J004859. 47— 733538.7 


TT~i A O T7 O O /(71 Tor 

IKAb 100471—735 


x-AGB 


—8.2 


Too f f\r\r\ to to *i 1 o r r 

JOUoOOU. 78 — 724125.5 




x-AGB 


—9.2 


J005030. 97— 730853.7 




x-AGB 


—9.7 


Tnorrino toocio et 

J 005100. 74—722518.5 


i\ /rev OA/T^ -\ CO 

MbA bMO lb3 


x-AGB 


—8.3 


Tnnc;i 1 o ct tqi nu^; /i 




A PD 
X-AbD 


—8.2 


J005446. 84— 731337.8 


T T7rT' i o r 

LlL/GC 105 


x-AGB 


—9.2 


Tnorrtroo no tootcto t 

J005530. 99 — 733752.7 




x-AGB 


—8.9 


Tnncrcr/io yi Tn,iTon o 

J UU5548. 54— 724729.0 




x-AGB 


—8.8 


TAnrrr /i f* o tomtom *i 

J005554. 63 — 731136.4 


KAW 960 


x-AGB 


—9.4 


JUuO ( lo. lz — f Z4ZOO.Z 




X-AbD 


— 0.3 


Tmri/ic;'? i q vori/iriQ o 
JU1U40O. lo — / zU4Uo.y 


OA A ACQ Tm fl/1 f^Q 


v APR 


— / .9 


JUU4oZ0.4o — / oZ04o.o 


OA A AQQ TfiH/l QO^J 


fA APR 


— 9.o 


JUU4yiZ. < — ( oZDZO.D 




fA APR 


1(1 ^ 
— 1U.O 


JUU4yo4.Ul — / oUoo / .o 




P APR 


— 9.o 


JUU4yoD.4D — ( oUoO ( .1 




(~a APR 




Tflfi/IQ^fl f>7 '70/100.1 Q 
JUU4yoD.U / — / 




P APR 


1(1 1 


JUUDU1D.4U — / ooUo4. / 




P APR 


1(1 1 


JUUD11D.41 — ( Z4Z00.U 




P APR 


1(1 (1 


TOO^I 07 Q7 VO/1/1/1G Q 
JUUD1Z / .o / — / Z444y.o 




P APR 


— 1U. z 


juudi^u.^o — / zo / zo.y 


A/T^IV tlMr 1 1 AO 
IVlO-A. OlvlV^ 1^±Z 


P APR 


Q ^ 


JUUODU ( .oo — ( olo4Z.O 




P APR 


1 (i 1 


J005617. 52— 722704.3 


O M TV !f f~A O o a o o o 

a°MO 204803 


C-AGB 


—9.8 


Tnncr ^ 1 no to^iooi n 

J 005o4 1.02 — 724831.0 




C-AGB 


—9.8 


Tnnc7Ci on to ,10,10 i 

J 005753 . 29 — 724343. 1 




C-AGB 


— 10.4 


Tnntrni o a-\ 7007 1 

J 0059 lo. 41 — 722734.2 




C-AGB 


—9.8 


Tnncnon nn too no n 

J005939. 00— 722308.0 




C-AGB 


— 10.3 


TO "1 00 OO O O 70007*77 /- 1 

J010029. 38— 722257.6 




C-AGB 


— 10.4 


Tnnoon 1 / ■ toooo ,1 "7 

J 003201. 63 — 732234.7 


T TT. r 1 1 OOO 

HV 11223 


O-AGB 


— 10.0 


Ton 1 /i ri 00 t n tj ct 1 1 ,1 

J 004249.83 — 725511.4 


T TT f 1 ocr 

HV 136b 


O-AGB 


— 10.1 


TOO A O A O 1 "1 TOO AAA O 

J 004843 . 1 1 — 730444. 




O-AGB 


— 10.7 


T O O •! O O "1 • > T T O O O O T O 

J004921. 37— 730327.2 




O-AGB 


— 10.0 


Tnndncn t/" 1 TO,iotrn 

J004950.7o — 724350.0 




O-AGB 


— 11.2 


Too/toro or 7QflQ1 ^ fl 

JUU4y52.U5— / 30316.0 




( \ a pri 


— 10.8 


TOOCOOT OO 70in01 O 

J005027. 20— 731021.9 




O-AGB 


— 10.8 


Tnntrooo 00 too 010 

J 005338. 83 — 733318.9 




O-AGB 


—9.2 


T r\nrz a a o c toocio ,1 

J005443.62 — 733512.4 




O-AGB 


— 10.1 


Tnnc cn 1 t Toioorr et 

J005850. 17—721835.5 


T TT T 1 Ol A O 

HV 12149 


O-AGB 


—9.0 


Tmnono a Tom co 

J 010302. 42 — 720153.0 




O-AGB 


—8.8 


TAI O ,1 O i • f TO O .4 ,1 O "1 

J010426. 65— 723440.1 


HV 1963 


O-AGB 


—9.8 


TOO^OO/1 ,1 O TO/IOMO /I 

J0049 24. 48 — 724312.4 




aO-AGB 


— 10.3 


Tnncnnn 01 toooeto 1 

J005000. 01 — 730853.1 




aO-AGB 


— 10.4 


TOO 7 00 7 OO TOOCOO 1 

J005005. 38 — 730500.1 




aO-AGB 


— 10.6 


TOOCOOT Ol TO ,10 00 o 

J005007. 31 — 724329.9 




aO-AGB 


— 10.4 


TOOP"OP"~l O"! TO O O O O O 

J005051. 81— 733036.8 




aO-AGB 


— 10.6 


Tnntri no £?o t o i i o ^ 

JUU51Uo.o2 — 724ol2.4 




aO-AGB 


— 10.4 


TOOr""IOO O A TO/IT'OT' /I 

J005109. 94— 724525.4 




aO-AGB 


— 10.5 


Tnnc ,in i ot toococt r 1 

J 00540 1.07— 733535. b 




aO-AGB 


— 10.7 


Toor"/ir"/i oo toooo/^ o 

J005454. 02— 730806.3 




aO-AGB 


— 10.4 


TOO 7 7J /I O OT TOOCTCTO ^ 

JUU5542.97— 73U552.D 




aO-AGB 


— 10.2 


JUU5o5o.oU— /2oyl4.o 




„ r\ A rA to 


— 10.2 


Toocrnc/i ncr toooi n r 1 

J 005954. 05 — 722219. D 




aO-AGB 


— 10.4 


Tnmo/ir o/^ tooooo t 

J00484D. 3d — 732820.7 




RSG 


—9.8 


TOO/IOCO TO TOnT/1£? o 

J004953. 78 — 730746.2 




RSG 


— 10.9 


J005006. 34-732811.0 


MSX SMC 096 


RSG 


-9.5 


J005021. 22-730609.5 




RSG 
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Fig. 2.— Sample SEDs for stars fit with GRAMS models. Pho- 
tometry from MCPS, 2MASS, and Spitzer is marked by filled 
squares. Near-IR photometry from the IRSF is marked by aster- 
isks, and AKARI photometry is marked by open triangles. Spitzer 
IRS spectra are plotted with solid lines, and the best-fit GRAMS 
model is marked with dotted lines connecting synthetic photometry 
at the observed wavelengths. 

stars may have gas-to-dust ratios similar to Gala ctic val- 
ues, even when in metal-poor en vironments (e.g.. lHabinel 
1996; iGroenewegen et al]|2007t ). For C-rich sources, we 
therefore assume ip = 200. 

We note that the GRAMS grid does not account for 
meta l lic iron dust, though recent work ([McDonal d et al.l 
120101 12011allbh shows that it may be an important con- 
tributor to dust production, especially for low-mass, 
metal-poor stars. The inclusion of metallic iron may 
provide substantial changes to the modeled DPRs for 
the fainter AGB stars; the direction of this change de- 
pends on competing factors relating to the higher opac- 
ity of metallic iron dust. Modeling of O-rich stars in 
the globular cluster u> Centauri suggests that the in- 
clusion of metallic iron in a purely dust-driven wind 
wi ll increase the DPRs by app roximately a factor of 
2 (jMcDonald et al.l 120091 1201 lbl ). However, we expect 
metallic iron to be only a minor component of optically 
thick dust shells, which dominate the dust production. 

2.3. Extrapolation to the Entire SMC 

In Figure [3l we show the DPRs as a function of 8- /mi 
excess for the stars fit in the previous section. The best 
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fits in Figure [3] are as follows: 
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Fig. 3. — 8-/im excess vs. D for stars marked in Fig. [l] DPRs 
were determined by applying the GRAMS model grid (see text). 
The solid lines mark the best fits for the filled symbols, and the 
dashed lines mark the best fits for the open symbols. In the middle 
panel, the dotted line marks the best fit obtained by combining 
the O-AGB and aO-AGB stars. We use the combined fit to derive 
the DPRs of the entire O-rich AGB sample to facilitate a direct 
comparison to the LMC a nalysis, which did not d istinguish between 
aO- and O-AGB sources (Srinivasan ct al. 200J3). 




3.0 3.5 4.0 4.5 5.0 5.5 6.0 3.0 3.5 4.0 4.5 

log[ Luminosity (L sun ) ] 

Fig. 4. — The DPR as function of luminosity for the SMC and 
LMC. The long-dashed line indicates the classical mass loss limit 
ll Jural 11984, D = L(t) cxp c)~ 1 <x L 75 ) scaled by the metallicity 
of the LMC (grey boxes) and SMC (black diamonds), the dash- 
dot line indicates the nuclear consumption rate for the CNO cycle, 
and the dotted line marks the nuclear consumption rate for the 3a 
process. Both nuclear consumption rates are directly proportional 
to L. LMC and SMC sources are plotted with light and dark points, 
respectively. Only sources with >3-cr excess are included. We note 
that RSGs with log[L/L0] < 4 are more likely to be AGB stars. 



relationship derived by Srinivasan et al.1 (|2009D . As with 
the SMC, the aO-AGB star DPRs are computed using 
the O-AGB relation, and the FIR object DPRs (assum- 
ing they are evolved stars) are computed using the x - 
AGB relation. The R SGs from Ivan Loon et all fl9991. 
Ivan Loon eFall (|2005fl , and iGroenewegen et alT(|2009bD 
are used to calibrate Equation [TJ resulting in: 



logl) (M yr- 1 ) = A + B log A 8 . Mm (mJy), (1) 
where 

A = -10.7, B = 1.2 (x-AGB/FIR), 
A = -10.4, B = 0.6 (C-AGB), 
A = -10.4, B = 1.1 (O-AGB / aO- AGB) , 
A = -10.7, B = 0.6 (RSG). 

The DPR of each evolved star in the SMC is estimated 
from its 8-/xm excess by applying the above Ag_ Mm -D re- 
lationship. Assuming the FIR objects are dusty evolved 
stars, their DPRs are derived using the x-AGB Xg-^m-D 
fit, as most of the FIR sources fall within the same mid- 
IR color-magnitude locus as the x-AGB stars. However, 
we stress that the majority of the FIR sources are likely 
Y SOs. 

ISrinivasan et al. (2009) considered only the x-AGB, C- 
AGB, and O-AGB candidates in their analysis of the 
LMC. Here, we recompute the DPRs for the LMC, sep- 
arating the aO-AGB and FIR sources and including the 
RSGs. These DPRs were computed using the X^.^-D 



A = -11.1, B = 1.2 (LMC RSGs). 

The resulting DPRs for both galaxies are plotted 
against luminosity in Figure 21 We see that all of the 
x-AGB stars lose mass at a rate that is higher than the 
nuclear consumption rate (dash-dot and dotted lines in 
Fig. [4j, implying that mass loss dominates their subse- 
quent evolution. The same is true for a subset of the 
C-AGB and O-AGB stars. 

2.4. The Cumulative Dust-Production Rates 

The cumulative DPRs are shown in Figure [5] and the 
total dust inputs are listed in Tables [2] and [3] The cu- 
mulative DPRs in both galaxies follow similar trends. 
The x-AGB stars dominate the total DPR at L > 4- 
8 x 10 3 Lq. The O-rich stars dominate the dust input at 
lower luminosities, and the regular C-AGB and the O- 
rich AGB stars (O-AGB + aO-AGB) contribute similar 
total DPRs (differing by a factor of 1 in the LMC and 1.6 
in the SMC). The FIR objects play a much larger role 
in the LMC, especially at luminosities higher than the 
classical AGB limit. This is likely due to increased con- 
tamination from YSOs in the LMC, which has a higher 
star formation rate. 
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Fig. 5. — Cumulative DPR vs. luminosity for the SMC (left) and LMC (right). The vertical dashed line marks the classical AGB limit. 



TABLE 2 

Dust-production rates in the SMC 



Type 


N 


ED 




% of 


< D > 








(M @ /yv) 




Total 


(M Q /yr 


) 


x-AGB 


313 


6.26 x 10" 


7 


65.9 


2.0 x 10" 


9 


C-AGB 


1559 


1.21 x 10" 


7 


12.7 


7.8 x 10" 


11 


O-AGB 


1851 


0.52 x 10" 


7 


5.5 


2.8 x 10" 


11 


aO-AGB 


1243 


0.26 x 10" 


7 


2.7 


2.1 x 10" 


11 


RSG 


2611 


0.31 x 10" 


7 


3.3 


1.2 x 10" 


11 


FIR 


50 


0.96 x 10" 


7 


10.1 


1.9 x 10" 


9 


Total 


7 627 


9.5 x 10" 


7 


100 


1.2 x 10" 


10 


Total, no FIR 


7577 


8.6 x 10" 


7 


90.5 


1.1 x 10" 


10 



due to an incorre ct assumption about the circumstellar 
dust composition. IGroenewe gcn (200(f) shows that if car- 
bon stars have less SiC and more amorphous carbon, for 
instance, the color-derived DPRs will decrease, bringing 
them closer to the excess-derived DPRs. SiC abundance 
is kno wn to decrease with metallicity (e.g. JLagadec et all 
so is expected to be uncommon in SMC stars. For 
O-rich stars, a lack of metallic iron or impurities in the 
silicates will decrease the 3.6-/im flux and a change in the 
shape of the 10-/zm silicate feature will affect the 8-/zm 
flux; both will affect how the [3.6] — [8] color scales with 
the DPR. Therefore, with the proper assumptions about 
dust composition, IR colors may provide DPRs similar 
to our excess-derived DPRs. 



TABLE 3 

DUST-PRODUCTION RATES IN THE LMC 



Type N ED % of < D > 
(M /yr) Total (M /yr) 

x-AGB 886 8.62 x 10~ 6 61.1 9.7 x 10 -9 

C-AGB 5190 0.87 x 10" 6 6.2 1.7 x 10" 10 

O-AGB 8871 0.63 x 10" 6 4.5 7.1 x 10" 11 

aO-AGB 6372 0.32 x 10~ 6 2.3 5.0 x 10" 11 

RSG 3908 0.24 x 10~ 6 1.7 6.1 x 10" 11 

FIR 179 3.40 x 10~ 6 24.1 1.9 x 10" 8 

Total 25 406 14.1 x 10" 6 100 5.5 x 10" 10 



Total, no FIR 25 389 10.7 x 10" 6 75.9 4.4 x 10" 



We note that stochastics play a role at high luminosity 
in Figure [5j since the brightest evolved stars are short- 
lived, and therefore rare. In the SMC, the 10 most ex- 
tr eme x-AGB s tars c ontribute 17% of the total DPR. 

iBover et al.l (|201 lh compute the dust input using the 
[3.6] — [8] color to compute the DPRs, assuming a dust 
composition of 85% amorphous carbon + 15% SiC for 
the carbon star s and 60% silicates + 40% A10 x for the 
O-rich stars (cf. Groenewegcn 2006). The overall trends 
in the total cumulative DPRs are the same, whether us- 
ing the IR excess or the IR color, but the color-derived 
rates are 4-7 times higher. This discrepancy might be 



2.5. Uncertainties in the Dust- Production Rates 

The results and analysis presented in this paper are 
very sensitive to the DPRs estimated by fitting the ob- 
served SEDs to GRAMS models. Some parameters, such 
as the outflow velocity, are unknown, so it is impossible 
to determine a formal uncertainty in the DPRs. In this 
section, we list the known uncertainties associated with 
the DPR m easurements . Some of th ese were already 
discus sed in I Sargent et al.l (|2011D and ISrinivasa n et al.l 
(|2011h . but we repeat them here for completeness. 

The current version of the GRAMS grid is computed 
for spherical dust shell geometry, assuming a constant 
DPR throughout the AGB lifetime. A treatment of the 
superwind phase is not considered at present. The DPR 
is computed from the dust shell inner radius, the size- 
and composition-averaged dust grain opacity and the 
outflo w velocity in the shell (cf. Eq. 2 inlSrinivasa n et all 
l2Ol0l ). When comparing the GRAMS-predicted DPRs 
with predictions from other studies, it is important to ac- 
count for the differing choices of dust optical constants; 
as mentioned in Section 12.21 

In GRAMS, the outflow velocity is assumed con- 
stant and i s fixed at 10 km s -1 , typical of metal- 
poor stars (Marshall et al. 2004). Here, we scale the 
DPR to account for the luminosity and the gas-to- 
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dust rati o (a nd hence metallicity) dependence of v eX p 
(Sec.lO e.g..lHabing et alJfl99i llyeiic fc Elitzurlll995l: 
van Loon et al.ll2008HWachter et al.ll2008l; lLagadec et al.l 
2010). 

In performing SED fitting to individual stars using the 
GRAMS models, it is necessary to take into account the 
uncertainties due to photometric errors as well as stel- 
lar variability. This results in a range of possible model 
fits, which can be translated to an uncertainty in the 
fit parameters, such as luminosity and DPR. A detailed 
SED-fitting study (D.Riebel et al., in preparation) finds 
that the relative error on the best-fit DPR is 10%-50% 
for the entire sample of significant mass-losing sources in 
the LMC. We expect a similar range of uncertainties in 
the DPRs used here, though this uncertainty is statisti- 
cal and diminishes by integrating over an ensemble. Also, 
the relative uncertainty in comparing the SMC with the 
LMC is not as large as the absolute uncertainty in each 
of their DPRs, owing to the similarity of the data (pho- 
tometry source, distance, and morphological type of the 
galaxies) and the treatment with the same model grid. 

The uncertainties discussed above relate to the deter- 
mination of DPRs for individual stars in our sample. The 
Xg-^m-D relations used in this paper are power-law fits 
to a representative sample of sources in the SMC (Fig. [3]), 
and there is some uncertainty in these fits. For instance, 
a slope that is 5% shallower for the Xg.^-D relation 
(such that B is 5% smaller in Eq. [1]) would result in 
DPRs 10-30% lower and result in a higher total contri- 
bution from stars with low DPRs. 

An additional source of uncertainty lies in our exclu- 
sion of binarity and the assumption of a spheroidal en- 
velope. In both cases, the inclination angle could affect 
the observed 8-/xm excess. It is impossible to determine 
binarity and/or geometry from IR photometry alone, but 
since this analysis includes a large sample of stars (pre- 
sumably with no preferred inclination angle with respect 
to the line-of-sight), this is not expected to be a large 
effect. 

3. DISCUSSION 
3.1. Circumstellar Dust at Low Metallicity 

In Tables [2] and we show the average DPRs for stars 
in the LMC and SMC, and in Figure IU we show how the 
DPR scales with luminosity. Some of the scatter in Fig- 
ure|4]is due to a spread in the evolutionary phase at each 
luminosity. LMC RSGs, which are exclusively O-rich, 
lose dust at a rate that is up to 10 x higher than SMC 
RSGs at the high-luminosity end (log[L/L Q ] > 4.5), con- 
sistent with the findings from several works showing that 
O-rich dust is more difficult to form at low metallicity 
(Ivan Loonl l2000l 120061: ISloan et al.l |2008). SMC RSGs 
with low luminosity appear to have higher DPRs than 
their LMC counterparts. This is likely due to uncertainty 
in the slope of the X^.^-D relationship (Sec. 12. 5[) . 

Previous studie s find that the DPR is lower in meta l- 
poor C-rich stars (jvan Lool][2000trvan Loon et al.ll2008h . 
On average, our findings agree with this, with mean 
DPRs 2.2-4.9 times higher in the LMC (Tables [2] and [3]). 
However, Figure 0] shows that the LMC C-rich stars can 
reach both high and low DPRs compared to the SMC, 
and it is clear that metal-poor C-rich AGB stars c an still 
be prolific dust producers (e.g.. ISloan et al.ll2009D . 



Using a simi lar treatment of the DPR, 

iGroenewegen et ail (|2007l ) find no clear metallicity 
dependence among Magellanic Cloud carbon stars. 
Since our individual DPRs are scaled by the luminosity 
and gas-to-dust ratio, they are a f actor of 3-9 lower 
than those from IGroenewegen et all (|2007l ) , though this 
should not change the relative DPRs in the SMC and 
LMC, since we use the solar gas-to-dust ratio (tp = 200) 
for the carbon stars in both galaxies. These results 
highlight the importance of AGB dust production at 
high redshift, where C-rich AGB stars can contribute 
strongly to the stellar dust production. 

3.2. Gas Input from Cool Evolved Stars 

The total current dust-injection rate from the cool 
evolved stars in the SMC is (8.6-9.5) x 10~ 7 Af Q yr -1 , 
depending on the number of FIR objects which are 
evolved stars. Assuming that the gas- to-dust ratio (iji) 
scales with metallic ity for all stars (|van Loonl [2000: 
iMarshall et "ail 120041) . then Vsmc = 1000 (see Sec. &2$. 
This yields a total gas return of (8.6-9.5) x 10~ 4 M Q yr -1 . 
Using Galactic metallicities for the C-rich objects (tp = 
200) and ip = 1000 for the O-rich stars (see Sec. [2T2|) 
yields a total gas return of (2.6-2.8) x 10 -4 M© yr -1 , 
which is close to the evolved star gas return estimates 
(derived from dusty evolved stars) in more metal-poor 
dwarf irregular galaxies with similar mass to the SMC 
such as WLM and IC 1613 (ZD = (3-7) x 10~ 4 M Q yr" 1 , 
or «10" 4 Ma yr" 1 if C-rich sou rces are adjusted to solar 
metallicity; iBover et al.ll2009bt ) . 

3.3. Other Dust Sources in the SMC 

We expect additional dust input from SNe and hot 
massive stars such as luminous blue variables (LBVs) 
and Wolf-Rayet (WR) stars, though the exact amount 
is highly uncertain. SNe can destroy dust as well 
as produce it. The SN dust destruction rate is also 
highly uncertain, and a large fraction of SN dust may 
not survive long enough to mix with the ISM (see 
Sec. I3.4j) . However, several SNe appear to harbor 
small amounts of dust. This includ es 0.02 M Pl of 
dust around SN2003gd i n NGC628 (|Sugerman et al.1 
[2001 lErcolano et al.ll2007h . 0.02-0.05 M Q around CasA 
(IRho et al.ll2008l), 2.2-3.4 x 10~ 3 M Q around SN2006jc 
in UGC4905 dSakon et al.1 [2009h. 1 xjC T 3 M Q in 
the Crab SN remnant (jTemim et all 120061 ). and 3 x 
10~ 3 M G) around the SM C SN remnant IE 0102.2-7219 
(jSandstrom et al.l [2009). Surprisingly, recent Herschel 
far-IR and ground-based sub-mm observations reveal a 
massive r eservoir of cold dust a r ound SN 1987A of 0.4 - 
0.7 M Q (|Matsuura et all 120111: lLakicevic etaLl 1201 ID . 
Based on these observations and excluding SN 1987A, 
we might expect SNe to produce ~ (0.1-5) x 10~ 2 M Q 
o f dust, on average. 

Filipovic et al.l (|1998t ) estimate the S N rate (tsn) in the 
SMC to be about 1 every 350 ± 70 yr. iMathewson et al.l 
(1983) estimate a more conservative rate of 1 SN every 
800 yr. These rates combined with the observed range 
of SN dust mass yields a SN dust input of (0.1-14) x 
10~ 5 Mq yr -1 , the lower limit of which is comparable to 
the input from AGB stars. 

The SMC is home to 3 LBVs and 5 supergiant 
B[e] stars that are detected at 24 am with MIPS 
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TABLE 4 
Dust Sources in the SMC 



Source 


Total Dust Input 




(M Q /n) 


C-rich AGB a 


(7.5-8.4) x 10" 7 


O-rich AGB 


7.8 x 10~ 8 


RSGs 


3.1 x 10" 8 


LBVs+B[e] stars b 


~io- 6 


SNe b ' c 


(0.13-14.3) x 10~ 5 



a The range in dust input from C-rich AGB stars 
depends on the fraction of FIR stars which are evolved 
stars. 

b The dust input from LBVs, B[e] stars, and SNe is 
highly uncertain. Much of the dust created by mas- 
sive stars may be ulti mately destroyed by the SNe 
shock. See Section [3731 

c The range in SNe DPR is computed assuming that 
SNe input ~ (0.1-5) x 10~ 2 M of dust, on average, 
and occur at a rate of 350-800 yr (Sec. I3.3[) . 

(jBonanos et al.ll2010D , which suggests the presence of cir- 
cumstellar dust. The dust production by B[e] stars is not 
well constrained. Dust inferred from the IR excess may 
be the remains of a debris disk left over from the forma- 
tion of the s tar rather than from dust formation in a stel- 
lar outflow. iKastner et al.1 (|2006f ) estimate the dust mass 
in the disks around two B[e] stars in the LMC, finding 
3 x 10~ 3 Af© of dust in the disk around R126. If formed by 
the star, this dust must be ejected in episodic/eruptive 
events since the mass-loss rate required is much too high 
to be sustained over the lifetime of a B[e] star. Each 
B[e] star might thus inject ~10 -3 Mq over its lifetime. 
The B[e] phase lasts ~10 5 yr, amounting to a total dust- 
injection rate of 10 -7 Mq yr -1 for 5-10 B[e] stars. This 
is similar to the dust input from C- and O-AGB stars in 
t he SMC. 

IBover et al.l ([2010D detect a dust-production rate of 2 x 
1O" 8 M yr" 1 around the LBV R71 in the LMC. If this 
rate is typical of dusty LBVs, then their contribution to 
the total dust budget is similar to the total DPR of B [c] 
stars. Dust production in LBVs is likely episodic, so if 
we assume that an LBV might lose a few solar masses 
of material over its lifetime (~10 4 yr) and assume a gas- 
to-dust ratio of 1000 for the SMC (see Section [372]), then 
each LBV might create 0.01 Mq of dust, corresponding 
to a total DPR of — 10~ 6 Mq yr" 1 and matching the 
input from the cool evolved stars. Galactic LBVs such 
as AG Carinae and r\ Carinae create 10 time s more dust 
than this due to a small er gas-to-dust ratio (|Voors et al.l 
120001: ISmithlTal1[2T)03l ). 

Table|4]summarizes the estimated dust input from each 
dust source. We note that much of the dust from massive 
stars may be destroyed in the ensuing shocks when these 
stars explode as SNe. Ultimately, the contribution to the 
total dust budget from massive stars may be negligible. 
We consider the full range of possibilities in Section 13.41 

3.4. Dust Survival in the ISM 

The dust input into the ISM from stellar sources is 
listed in Tablegl The dust mass in the SMC ISM inferred 
from IR and sub -mm imaging is (0.29-1.1) x 10 6 Mq 
(|Bot et all [2010h . The likelihood that the ISM dust 
is stellar in origin depends on the lifetime of dust in 



the ISM. Several works find that SN shocks are the 
dominant dust-destruction mechanism in t he ISM (e.g., 
iDraine fc Salpeterl Il979t Uones et al.l I1994D . More than 
50% of a silicate grain (a = 0.1 ^m) is returned to the gas 
phase i n shocks with v s > 200 km s _1 (|Draine fc Salpeterl 
1979). lDT7ah^(l2009l) argue that a typical SN with energy 
Eq = 10 51 erg in a medium with density njj = 1 cm~ 3 
will remain in the Sedov- Taylor phase until the shock 
speed drops to 200 km s _1 , so we expect a SN to process 
M w E /v% = 1260 M Q of interstellar material. The 
dust lifetime within the ISM is then: 

M gas 

= ISM /n\ 

d (1260 M /r SN ) ' 1 ' 

where M?™„ (H I + H 2 ) is 4.5 x 10 8 M (jBolatto et al.l 
[2?lll . I.Tones et all (Jl996) argue that if the grains are 
porous, their lifetimes could be enhanced by a factor of 
3. This , com bined with Equation [2] and using t§n from 
Section 13.31 results in an expected dust lifetime of Td = 
0.38-0.86 Gyr, similar to the lifetime estimate d for the 
Milky Way (-0.5 Gvr: Uones et~aT][T99l [1991 . 

It is possible that many SNe will explode in much 
denser molecular clouds, with nn — 10 3 cm~ 3 . In this 
case, the Sedov- Taylor phase terminates at higher shock 
velocities, and the SN may only process ~170 M of 
material, leaving much of the ISM dust unaffected. If all 
SNe exploded in dense molecular clouds, then the dust 
lifetime increases to Td = 2.8-6.7 Gyr. 

If no stellar dust is destroyed by SNe (or other mech- 
anisms such as grain-grain collisions), then the life- 
time of dust in the ISM is determined by the SFR as 
Td = A/ism/ SFR, wh ere SFR is the star formation rate 
(3.7 x 10" 2 Mq vr-^ lBolatto et alJ[20Tl . This yields a 
lifetime of r™ ax = 12 Gyr, very near the age of t he oldest 
stellar population in the SMC ([Noel et alJl2009| ). A dust 
lifetime of 12 Gyr would suggest that most dust that has 
ever been created in the SMC survives in the ISM today. 

In order to compare the dust input to the observed 
ISM dust mass, we now consider three cases: 

(A) massive stars ultimately contribute no dust (it is all 
destroyed by the SN shock), such that cool evolved 
stars are the only stellar dust sources, 

(B) progenitor dust survives the SN shock and SNe 
themselves produce the lower limit of the esti- 
mated DPR listed in Table H and 

(C) progenitor dust survives the SN shock and SNe 
themselves produce the upper limit of the esti- 
mated DPR listed in Table g] 

Based on the total dust input and the observed ISM dust 
mass, we can estimate the required dust lifetime for each 
scenario. 

For case A, the total dust input is (8.6-9.5) x 
10~ 7 Mq yr -1 and the dust must therefore survive for 
(3-13) x 10 11 yr, depending on tsn, to account for the 
observed ISM dust mass. This is 2-3 orders of magnitude 
longer than the lifetime estimated with Equation^ and 
indeed is longer than a Hubble time. Case A is therefore 
an impossible scenario unless the SMC underwent peri- 
ods of much higher star formation prior to the era that 
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created the current popula tion of AGB stars . In 12 re- 
gions spanning the galaxy, No el et all (|2009l ) show that 
the SFR in the SMC remained constant within a factor of 
3 from the onset of star formation to intermediate ages, 
so we do not expect a significantly larger contribution 
from AGB stars in the past than is currently observed. 

For case B, the total dust input is (3.3-3.4) x 
10 -6 Mq yr" 1 , so the dust lifetime must be (9-33) x 
10 11 yr. This lifetime is also much too long, so the ISM 
dust cannot be supplied solely by AGB stars and SNe if 
SNe produce only 10" 3 M , on average. 

For case C, the total dust input is 1.4 x 10 4 M Q yr -1 
and the dust lifetime must be 2.1-7.9 Gyr, comparable 
to the lifetime of porous grains if SNe Shockwaves are 
the dominant destruction mechanism in the SMC and 
tend to explode within dense molecular clouds (Eq. [5J 
r d = 2.8-6.7 Gyr). SNe that produce 5 x 1(T 2 M Q and 
occur every 350 yr can thus feasibly produce the observed 
ISM dust mass. We also note that if SNe can produce 
the amount of dust recently observed around SN 1987A 
(0.4-0.7 Mq), they can explain the ISM dust even for 
the most pessimistic estimates of the dust lifetime. 

3.5. The SMC Dust Budget: Excess Dust in the ISM? 

It is clear that AGB and RSG stars alone cannot ac- 
count for the observed ISM dust, even if SNe shocks de- 
stroy no dust at all. With certain assumptions, SNe and 
cool evolved stars together can be made to account for 
the observed ISM dust mass. This includes assuming 
that a high fraction of dust created by SNe and their 
progenitors survives and that most SNe explode in a 
dense medium, such that most of the surrounding ISM 
dust is shielded from destruction. These assumptions are 
somewhat generous, and we are left with at least some 
excess of ISM dust, si milar to what is seen in the LMC 
(jMatsuura et alll2009ft . the Milky Way (|Dwekll 19981 ). and 
in high-redshift galaxies. Assuming the dust lifetime de- 
rived from Equation[2l the AGB and RSG stars combined 
can account for only up to 2.1% of the ISM dust. 

Case B in Section I3T41 assumes a conservative estimate 
for the SNe dust production. This scenario impl ies that 
dust m ust grow in the ISM itself, as discussed by iDraind 
(2009) and references therein, unless evolved stars pro- 
duce more dust than is implied by the mid-IR obser- 
vations. Far-IR imaging of a central 2° x 8° strip of 
the LMC with Herschel revealed a strong far-IR ex- 
cess (> 200 jim) around only one star, the LBV R71 
(|Bover et a l. 2010), suggesting that cold dust envelopes 
are rare around evolved stars. Cool evolved stars such as 
AGB and RSG stars are thus unlikely to be the solution 
to the missing dust problem. However, if SN 1987A is an 
anomaly or destroys its own dust, the dust input from 
cool evolved stars, especially extreme carbon-rich AGB 
stars, rivals that from SNe. 

Dust accretion in molecular clouds migh t contribute 
significantly to the to tal dust budget fe.g.. lDweklll998l ; 
IZhukovska et al.l 120081 ). A chemical evolution mo del of 
the Milky Way derived by IZhukovska et all (|2008l ) finds 



that the AGB stars dominate the dust input until the 
metallicity surpasses Z 10~ 3 . The SM C metallicity 
is Z = (2 ± 0.7) x 10" 3 (jLuck et al.lll99l . so we might 
expect the dust input from AGB stars to rival the rate 
of dust growth in molecular clouds. Assuming case B in 
Section [3.41 we find that dust grains must grow at a rate 
of (0.1-3) x 10~ 3 Mq yr" 1 , using conservative estimates 
for the dust lifetime (0.4-2.2 Gyr). This is similar to 
the rate of dust growth in t he Milky Way ISM suggested 
bv lJenniskens et all (|1993l ) (10" 3 M yr" 1 ). The rate 
of dust accretion in the ISM may thus be 2-3 orders of 
magnitude larger than the AGB dust production. 

Alternatively, we must also consider that the ISM recy- 
cling timescale may be long enough such that the current 
stellar dust input is not representative of the stellar input 
during the epoch that created today's ISM dust. During 
an increase in star formation, the SNe and RSGs w ould 
contribute more to the dust input. iNoel et all ([2009D find 
that the SFR was enhanced 0.2-0.5 Gyr ago in the east- 
ern fields, perhaps corresponding to a close encounter 
with the LMC, and resulting in a larger rate of SNe dust 
injection during that era. 

It is also possible that the ISM dust mass is overes- 
timated. However, the dust mass assumed here ((0.29- 
1.1) x 10 6 M : IBot et"aTll2010D includes the sub-mm ex- 
cess, which is well-constrained by invoking spinning dust 
grains. If, instead, the sub-mm excess is due to very cold 
dust grains, the ISM dust mass would increase dramati- 
cally, exasperating the problem. 

4. SUMMARY OF CONCLUSIONS 

We estimate the dust-production rates of cool evolved 
stars in the SMC using the 8-/im excess and compare the 
dust input to that from SNe and to the dust mass in the 
ISM. We find that the C-rich AGB candidates account 
for 87%-89% of the total cool evolved star dust input. 
The equivalent fraction in the LMC is 89%-91%. The 
majority of this dust input comes from the extreme AGB 
stars. RSG stars play a minor role in the dust input, 
especially below the classical AGB luminosity limit. 

While we can now quantify the dust production in the 
winds of cool evolved stars, the SNe dust-production rate 
(and dust-destruction rate) remains poorly constrained. 
It is possible that SNe can account for all of the ISM 
dust if they can each produce the upper range of dust 
masses observed around SN remnants. If, on the other 
hand, SNe can produce only the smaller dust masses in- 
ferred from mid-IR observations of several SN remnants, 
we expect SNe to contribute equally to the ISM dust 
compared to the cool evolved stars. This is similar to 
the findings in the LMC, suggesting only a small vari- 
ance with metallicity. In this case, an additional dust 
source is required and perhaps implies that dust grows 
efficiently in the ISM. 
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